Prepared for submission to JHEP 



Hadronically decaying color-adjoint scalars at the LHC 



Steffen Schumann,"'^ Adrien Renaud,*^ Dirk Zerwas'^ 

Institut fiir Theoretische Physik, Universitdt Heidelberg, 69120 Heidelberg, Germany 
^ n. Physikalisches Institut, Universitdt Gottingen, 37077 Gottingen, Germany 
"LAL, IN2P3/GNRS, Orsay, France 

E-mail: s . schumannSthphys . uni-heidelberg . de, renaudOlal . in2p3 . f r, 
zerwasSlal . in2p3 . f r 

Abstract: We study the phenomenology of the pair-production of scalar color-octet elec- 
troweak singlet states at the LHC. Such states appear in many extensions of the Standard 
Model. They can be pair-produced copiously at the LHC and will signal themselves as 
resonances in multijet final states. Beyond the QCD pair-production process we consider 
a vectorlike confinement scenario with an additional color-octet vector state. These vector 
particles can be produced in the s-channel and through their decay contribute to the scalar 
pair production. We point out the differences between the two hypotheses and device a 
strategy to distinguish them. 
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1 Introduction 



The potential to produce new colored particles is one of the most prominent features 
of hadron colliders such as the Fermilab Tevatron or the CERN LHC. There, new-physics 
states that directly couple to the partonic initial state can be produced copiously for a wide 
range of masses. Most candidate theories for physics beyond the Standard Model (SM) 
predict the occurrence of such new particles and there is a continuous effort to search for 
hints for their production. Recent examples include the search for supersymmetric squarks 
and gluinos [1-4] or lepto-quarks [5, 6]. The various searches might be classified according 
to the considered final states. Besides a certain number of hadronic jets, remnants of the 
particle's color charge, this includes identified leptons from electroweak decays, associated 
photons and/or missing transverse energy from undetected decay products, e.g. a potential 
dark matter particle. 

In this paper we analyse the pair-production process for new scalar color-octet states 
that are singlets under the electroweak gauge group. Such states appear in various exten- 
sions of the SM either as new fundamental particles [7-9] or composite objects [10, 11]. As 
we consider states in the color-octet representation their QCD production cross section is 
fixed. The pair-production process, however, might have contributions depending on the 
theory that embeds the new scalars such as additional new particles that are produced 
resonantly in the s-channel and can decay into a pair of such octet scalars. The decays of 
the new scalar particles are of course highly model dependent and closely related to the 
underlying theory's spectrum. However, most likely, for a certain mass range two-body 
decays to SM partons will dominate and it is this decay mode we focus on here. As a 
consequence we have to search for the states under question in multijet final states, where 
they would signal themselves as dijet resonances. 

The precise signature for the hadronically decaying resonances certainly depends on 
the mass of the new states as well as the experimental trigger requirements. The con- 
tinuum QCD production process will favour pair production at the kinematic threshold. 
In consequence we consider inclusive four-jet production as the search environment. In 
situations where the octets get produced with sufficiently large transverse momentum a 
dedicated subjet analysis in dijet final states might be the most powerful search strategy 
[12]. 

The central challenge of the analysis is to dig up the signal from the enormous QCD 
multijet background that exceeds the signal by orders of magnitude. The key point is to 
make optimal use of the kinematics of the signal events to represent the production of two 
equal mass narrow scalar resonances. In an earlier study the feasibility of this analysis 
was shown for the vector like confinement scenario [11]. Here we further improve that 
analysis, extend it to discovering the pure QCD pair-production process, and point out how 
different models could be differentiated experimentally. We propose a data-driven approach 
to extract the relevant QCD multijet background, thereby minimizing the dependency on 
theoretical predictions. 

In principle our candidate particles can also be produced singly at hadron colliders 
and there exist very stringent bounds from dijet resonance searches. For example ref. [13] 
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quotes a 95% CL exclusion limit of 1.91 TeV for a color-octet scalar resonance with an 
order unity coupling to two gluons [14]. However, as soon as this interaction is induced at 
the loop-level only, the single-production rate drops significantly thus evading the strong 
limit. Such scenarios we want to consider here, for which the pair-production process is 
the most promising discovery channel. 

The paper is organized as follows. In Section 2 we define our two benchmark scenarios, 
namely the most simple vectorlike confinement model and the SM extension by a complex 
scalar gluon that decays into two gluons with probability one. In Section 3 we describe 
the tools used to simulate signal and background events. In Section 4 we give details 
on our event-selection criteria and discuss a method that allows an extraction of the QCD 
background from data. In Section 5 we present results for our proposed search showing how 
one could potentially distinguish between the two benchmark models. Section 6 contains 
our conclusions. 

2 Two benchmark models 

Throughout this paper we want to consider two benchmark models for the production 
of electroweak singlet color-octet states. The first one is a vectorlike confinement scenario 
featuring besides our scalar candidate, the real octet hyper-pion, a further vector resonance 
in the color adjoint, called hyper-rho or coloron, that predominantly decays into a pair of 
hyper-pions. As the second template we consider the scalar partners of the gluino (called 
sgluons) as they appear in extended versions of the Minimal Supersymmetric Standard 
Model, such as the Minimal i?-symmetric Supersymmetric Standard Model (MRSSM) [15] 
oraniV = l/iV = 2 hybrid model of supersymmetry with the N = 2 supersymmetry 
implemented only in the gaugino sector [9]. In the following we will briefly deflne our two 
model hypotheses and discuss the different predictions for the pair-production rate of the 
two scalar octet candidates, the hyper-pion and the sgluon respectively. 

2.1 Hyper-pions in vectorlike confinement gauge theories 

The approach of vectorlike confining gauge theories at the TeV scale, first proposed in 
[16] and later on generalized in [17, 18], provides an attractive ansatz for potential physics 
beyond the SM. Assuming a minimal extension of the SM by new matter fields and making 
use of a fundamental mechanism already realized in nature, namely a confining strong 
force, a very rich phenomenology can be obtained from those theories described through 
rather few parameters only. The starting point is a quite minimal extension of the SM by 
new fermions in vectorlike representations of the SM gauge groups. These new fermions, 
charged under the strong and/or electroweak gauge group, feel a new strong gauge force 
called hypercolor (HC) that confines at the TeV scale, just as QCD confines at scales 0(100 
MeV). The emerging hyper-hadrons can be produced at the LHC signalling themselves 
through their decays mediated by hypercolor or SM interactions. The LHC phenomenology 
of vectorlike confinement will be dominated by the production of the lightest scalar and 
vector states, namely the pseudo-Nambu-Goldstone bosons called hyper-pions n and the 
counterpart of the QCD p meson, the hyper-rho p. 
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Figure 1. Feynman diagrams for resonant coloron production decaying into a pair of color- adjoint 
scalars. Note that we consider e = 0.2 while we set ^ = such that the gluon-gluon fusion 
contribution vanishes. 

We want to examine the most minimal vectorlike-confinement model that has been 
studied before in refs. [11, 16, 19, 20]. We consider three new massless vectorlike fermions 
charged under QCD and hypercolor only. The hypercolor gauge group is fixed to SU(3)hc 
and we assign the fermions as bi-fundamentals of the SU(3)c (8'SU(3)hc product group and 
singlets under the electroweak interactions. The lowest lying states below the hypercolor- 
confinement scale Khc ~ TeV are the color-adjoint pseudo-Nambu-Goldstone hyper-pions 
vr and the octet hyper-rho or coloron p. For this scaled-up version of QCD their mass ratio 
is approximated by 



with rrip ~ Ahc- The effective Lagrangian describing the dynamics and interactions of the 
hyper-pion and hyper-rho is given by [11] 



The first line contains the SM and hyper-rho kinetic terms and the coupling of the 
hyper-rho to quarks induced by kinetic mixing of the gluon and coloron fields. The second 
line summarises the kinetic and mass terms for the hyper-pion and the effective vertices 
that describe the decays p — )• tttt and tt — >• gg. The terms in the third line represent strong 
interaction matrix elements of the underlying theory that cannot be extracted from the 
QCD analog [11]. The parameter ^ is an undetermined number of order one. We use the 
most conservative assumption of ^ = as nonzero values enhance the pion pair-production 
cross section, cf. ref. [11], through an additional resonant coloron contribution in the gluon- 
gluon fusion channel, cf. Figure 1. The term proportional to x is of no relevance for our 
analysis as it contributes to the coloron pair-production process only [11]. 

All other parameters of the model are fixed by just scaling up the QCD spectrum and 
the corresponding couplings of the low-energy QCD chiral Lagrangian [16]. The mixing 
parameter e is determined to e ~ 0.2 and the coloron to hyper-pion coupling is given by 
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Figure 2. Feynman diagrams for the QCD pair production of color adjoints in the qq and gg 
channel. 

g'p^jf- ~ 6. The hyper-pion decay constant is determined through 



Ahc Aqcd 



(2.3) 



accordingly the hyper-pion decays promptly into two gluons. The coloron is rather broad 
state, 

^«0.19. (2.4) 

Its dominant decay mode however is not into dijets but rather a pair of hyper-pions, 
determined by 

^P^m e (2.5) 



Tp— >7r7r 9pTTTT 



2.2 Sgluons in supersymmetric models with Dirac gauginos 

Supersymmetric models that accomplish Dirac gauginos provide an attractive solution to 
the flavour problem of the MSSM [15]. This is achieved through an additional suppression of 
dangerous flavour- violating processes. Considering a Dirac gluino, in order to account for its 
four degrees of freedom, besides the gluon two additional bosonic degrees of freedom must 
be added in a viable supersymmetric theory. This complex color-adjoint scalar, known as 
the sgluon, has explicit realisations in the context of the MRSSM [8] or a phenomenological 
N = \/N = 2 SUSY hybrid model as considered in [9]. 

The sgluon has a gauge coupling to gluons and its pair-production rate is thus given 
by QCD. It interacts with quarks only at the one-loop level with the couplings being 
proportional to the quark masses. The sgluon-gluon-gluon interaction also occurs at the 
one-loop level only and is mediated by squarks. Accordingly resonant single-sgluon pro- 
duction is strongly suppressed. The sgluon coupling to two gluinos is the supersymmetric 
partner of the gluon coupling, while the coupling to squarks originates from D terms and 
is proportional to the Dirac gluino mass [8]. 
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Figure 3. (Left) The tree-level total cross sections at the 7 TeV LHC for the production of a pair 
of octet scalers. Besides the hyper-pion signal we present the cross section for producing a pair of 
scalar gluons [8, 9]. (Right) The ratio of the hyper-pion and sgluon cross section as function of the 
scalar mass. 

The best laboratory to discover the sgluon is its pair-production process. In [8] two 
split sgluon states have been considered. The work focused on the spectacular decay of 
the sgluon into top or anti-top accompanied by a light-flavour quark what leads to a jet- 
associated same-sign top signature. If heavy enough the sgluon can decay at the tree-level 
into pairs of gluinos or squarks. It is those decays that ref. [9] focused on along with the 
decay into ti, however, the authors considered a degenerate pair of sgluons. 

In both scenarios the sgluon has a sizeable branching fraction into gluons that dom- 
inates certainly below the top threshold. It is this part of the sgluon phenomenology we 
want to study here. To be precise, we consider two degenerate sgluon states as in ref. [9] 
and fix the sgluon branching ratio to gluons to unity independent of its actual mass. Such 
we simplify the underlying extended supersymmetric model but we still capture a large 
part of its phenomenology. In order to establish the supersymmetric theory that embeds 
the scalar adjoint the decay channels discussed in refs. [8, 9] need to be studied. 

2.3 Octet-scalar pair production 

In Figure 2 we depict the Feynman diagrams for the QCD pair production of hyper-pions 
and sgluons from qq and gg initial states. The way we set up our sgluon benchmark 
model its pair-production rate exceeds the continuum QCD production rate for a pair of 
hyper-pions by a factor of two, consequence of the degeneracy of the two sgluon states. 
However, the actual hyper-pion pair-production rate is determined by the additional s- 
channel coloron contributions depicted in Figure 1. In Figure 3 (Left) we compare the 
7 TeV LHC cross sections for the hyper-pion and the sgluon pair-production signal as a 
function of the scalar mass. The ratio of the two cross sections is shown in Figure 3 (Right) 
as function of the mass. For low scalar masses the ratio is almost 0.6 as expected for a real 
scalar with respect to a complex scalar field. At higher masses the ratio approaches 0.83 
as the contribution of the s-channel coloron exchange becomes more important before 
dropping back to 0.6 at 1 TeV. 

The two models outlined above shall serve as concrete realisations of new color-adjoint 
scalars. As mentioned before such particles appear in many Standard Model extensions. 
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The pair production of electroweak singlets decaying into heavy quarks has for example 
been considered in refs. [21-23], rare decays have been studied in ref. [24], phenomenological 
analyses of electroweak-charged adjoint scalar production can be found in refs. [25-27]. 

3 Simulation of signal and background events 

Before turning into the discussion of our search strategy and signal selection criteria we 
want to briefly describe the simulation tools employed to model the signal and background 
processes. 

We use the Sherpa event generator [28] to simulate both, the hyper-pion and sgluon 
pair-production signal as well as the QCD multijet background. The signal models have 
been implemented in the matrix-element generator AMEGIC++ [29]. The Feynman rules for 
the hyper-color model have been derived from Eq. (2.2). The interactions relevant for the 
pair-production of sgluons are given by QCD and we assume the sgluon branching ratio to 
gluons equal to unity. 

For the signal events we set the factorisation and renormalisation scales equal to the 
scalars mass, i.e. fip = fJ-R = "^Ig- The decay of the scalar resonances into gluons we treat 
in the narrow-width-approximation, whereas we fully take into account off-shell effects 
for the broad coloron s-channel contribution in the hyper-color scenario. The parton-level 
signal events then get passed to Sherpa's parton-shower and hadronisation routines [28, 30] 
in order to obtain realistic particle-level events. 

The modelling of the QCD multijet backgrounds is a severe challenge [31]. Up to now 
the next-to-leading order corrections to four-jet production are unknown so we have to rely 
on leading-order predictions. To accurately simulate multijet production we employ matrix- 
element parton-shower matching as implemented in Sherpa [32, 33]. We consider the 
complete sets of tree-level matrix elements with up to six final-state partons from Sherpa's 
matrix-element generator COMIX [34]. The matrix elements of varying multiplicity get 
consistently combined to an inclusive sample of fully exclusive events through subsequent 
parton-shower evolution [30] before subject to hadronisation. The renormalisation and 
factorisation scales are dynamically determined on an event-by-event basis according to the 
matching algorithm [32] . For the parton-separation parameter of the matching procedure 
we use Qcut = 30GeV. Note that Sherpa's model for QCD multijet production has already 
successfully been validated against data from HERA [35], Tevatron [36, 37] and the LHC 
[38, 39], what makes us confident to use a realistic background estimate throughout this 
study. In Sec. 4 we describe a strategy to extract the normalisation and shape of the QCD 
background in a data-driven approach. 

To account for detector-resolution effects such as finite granularity and jet smearing we 
process the particle-level signal- and background events through the publically available 
detector-simulation package Delphes [40]. Delphes yields a reliable simulation for a 
prototypical LHC detector that has been used in several physics studies such as [41-43]. 
For defining jets we make use of the FASTJet [44] package. 
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Figure 4. Distribution of the reconstructed scattering angle in the center-of-mass frame of the 
reconstructed four highcst-p^ jets. Note, for this shape comparison the signal was scaled by a 
factor of 10. 

4 Event selection and background estimation 

Throughout the analysis we use anti-fcy jets [45] with an R parameter of 0.6. Events are 
required to have at least four jets with a minimal px of 50% of the mass of the resonance. 
Once this pT selection is applied the signal is composed of octet scalars produced with a 
non-negligible transverse boost. In the i] — (j) plane the separation of the gluon jets from 
the decay is typically AR ~ 1. We then consider all the possible pairings between the 
four highest-pr jets and retain the combination that minimizes \ARij — 1| + {ARi^i — 1| 
where k, I denote the four leading jets. If one of the reconstructed resonances has a jet 
separation larger than ARij = 1.6, the event is rejected. 

To improve further the rejection of the background, the relative difference between the 
two reconstructed masses is required to agree within 7.5%, i.e. \Mi — M2\/{Mi + M2) < 
0.075. Additional separation between background and signal is obtained by using the 
reconstructed scattering angle (cos(0*)) in the center-of-mass frame of the four jets to be 
less than 0.6. The cos{9*) distribution is shown in Figure 4 for the hyper-pion signal and 
the QCD background. All cuts but the cut on the scattering angle were applied. For 
the signal the distribution is sin(2^*) like as the production of scalar particles is central, 
whereas for the QCD background the production in the forward region is more pronounced. 

The mean of the two reconstructed masses for events that pass all cuts is shown in 
Figure 5 for a hyper-pion signal of m^p = 100 GeV. The data corresponds to an integrated 
luminosity of lfb~^ at the 7 TeV LHC. As the form of the signal is similar to the form 
of the background, the ABCD matrix method is used to predict the background. The 
four regions are separated via the scattering angle (less than or greater than 0.6) and the 
relative mass difference (inverting the cut described above), cf. Figure 6. The region with 
the highest signal purity for that we want to estimate the background normalisation and 
shape is denoted as A. Under the assumption our two discriminators are uncorrelated for 
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Figure 5. The average mass for a hyper-pion signal of ma- = 100 GeV, the QCD background and 
the sum of signal plus background are shown after all cuts for an integrated luminosity of lfb~^. 
The green line is the predicted background according to the ABCD method. 
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Figure 6. The four event categories used in the ABCD method. The two discriminators are 
the scattering angle cos(6'*) and the relative mass difference of the two reconstructed resonance 
candidates. The signal-enriched region corresponds to selection A. 



the QCD background, we can estimate the number of background events in region A to 

(4.1) 



BG 



The background shape we take from region D, where only the cut on the scattering angle is 
reversed with respect to the signal selection. The background prediction from the ABCD 
method is shown in Figure 5. Both the background shape as well as the normalisation 
are predicted well when comparing to the Monte-Carlo truth. This data driven approach 
will have to be validated with actual collision data of course, i.e., will the same correlation 
pattern be observed in the data as predicted by Monte Carlo and which scale factor will 
have to be applied to the Monte-Carlo prediction to account for higher-order corrections. 
Here our emphasis was on the background per se, in the actual data potential signal 
contaminations in regions B, C and D will have to be taken into account as well. 
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Figure 7. Cross sections after all cuts for the QCD background, and the hyper-pion signal. 
5 Results 

In Figure 7 the cross sections after all cuts are shown as function of the scalar mass. The 
QCD background decreases strongly as expected. 

In Figure 8 the sensitivity of the analysis is shown as function of the hyper-pion mass. 
We define the sensitivity as the statistical error on the background (and its prediction) 
for lfb~^ in a window of 40 GeV centered on the nominal hyper-pion mass divided by 
the efficiency and multiplied by 1.64 (similar to a one-sided Gaussian counting experiment 
limit). 

As the instantaneous luminosity increases the trigger thresholds of the experiments 
will increase, using as we are using a sliding pT cut, cf. Sec. 4, the analysis is valid for 
masses greater than twice the trigger threshold. Lower masses than that are not accessible. 
For masses far from the trigger threshold one could perform a sideband analysis only in the 
signal region (bump hunting), i.e., an analysis with a fixed cut on the transverse momentum. 
This however can lead to double peak structures for the signal, diluting somewhat the gain 
on the efficiency. 

For a mass of 300 GeV the hyper-pion (sgluon) cross section is 10 pb (13 pb). The 
sensitivity of the analysis is 5 pb. Thus even adding experimental errors, the unknown 
K-factor of the background cross section, we can surely state that the four-jet final state 
can be probed in spite of the huge QCD background. 

The difference between the hyper-pion and sgluon pair production goes beyond a scal- 
ing of the cross sections. The s-channel coloron exchange also has an impact on the 
kinematics and therefore on the selection efficiency. In Figure 9 the transverse-momentum 
distribution of the 4th jet is shown for a signal mass of 300 GeV. The s-channel coloron 
increases the pT and therefore for the same analysis the efficiency for hyper-pions and sglu- 
ons is not the same as shown in Table 1. While the efficiency for the hyper-pions is stable, 
the efficiency for sgluons decreases slowly for increasing masses. It is interesting to note, 
that the efficiency as a function of the pT cut behaves differently in the two scenarios. For 
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Figure 8. The sensitivity of the analysis is shown, based on the statistical error in a window of 
40 GeV centered on the nominal mass. 
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Figure 9. The distribution of the 4**^ jet is shown for the sgluon and hyper-pion signal for a scalar 
mass of 300 GeV (before cuts). 

example for a scalar mass of 500 GeV a decrease of the pT cut by 20% leads to an increase 
of the hyper-pion selection efficiency by a factor of 2, the efficiency of the sgluon signal, 
however, increases by a factor of 3. A dedicated additional analysis for masses far from the 
trigger threshold could be envisaged that could provide further information to discriminate 
the two models. 

100 GeV 225 GeV 300 GeV 500 GeV 



hyper-pion 0.9% 1.1% 1.2% 1.2% 
sgluon 0.6% 0.4% 0.3% 0.2% 



Table 1. Efficiencies are shown as function of the scalar mass. 
An interesting question arises also whether the coloron mass can be reconstructed. To 
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Figure 10. The distribution of tlie reconstructed coloron mass is shown after all cuts for the 
hyper-pion model and the sgluon model. For the latter the distribution is flatter, whereas for the 
former the s-channel contribution is visible. From left to right the distributions correspond to scalar 
masses of 100, 300 and 500 GeV. In the hyper-color scenario the corresponding coloron masses are 
333, 1000 and 1667 GcV, respectively. 



address this question we consider the invariant mass of the two reconstructed scalar can- 
didates. As can be seen in Figure 10 at low mass the coloron mass distribution is diluted 
by the non-resonant production mechanism. No clear distinction can be observed between 
the sgluon and hyper-pion case. However for masses of 300 GeV as well as 500 GeV, cf. 
the middle and right panel of Figure 10, the coloron of the hyper-pion model becomes 
visible, thus giving us possibly a handle to distinguish the two models. We note however 
that in order to differentiate the two hypotheses, the actual QCD background (larger than 
the signal) will have to be dealt with. While the signal differences are quite striking, an 
integrated luminosity much larger than lfb~^ will be necessary to conclude on the under- 
lying model and only if the signal turns out to be at sufficiently high mass. However, the 
simulation of such large integrated luminosities for the background is beyond the scope of 
our paper. Furthermore we want to add that for an LHC center-of-mass energy of 14 TeV 
the discrimination of the two models through the reconstruction of the coloron resonance 
gets more complicated due to the different initial-state parton fluxes. For the hyper-pion 
mass range considered here, namely = 100. ..500 GeV, at 14 TeV the hyper-pion to 
sgluon cross-section ratio is strictly smaller than for 7 TeV. In turn the relative importance 
of the coloron contribution to the total hyper-pion pair-production cross section is largely 
reduced and the QCD continuum production mechanism dominates. On the other hand 
this dependence on the collision energy might serve as a handle to establish the coloron 
production from qq initial states. 

The alternative way to verify the existence of the hyper-rho state is to study its pair- 
production process. In refs. [11, 19, 20] the signature — )• — )■ 47f — )• 8jets has been 
studied and it was shown that the LHC provides great discovery potential both at 7 and 
14 TeV center-of-mass energy. 



6 Conclusions 



We have presented a search strategy for the pair production of new color- adjoint scalar 
particles that decay into pairs of jets. Such states are predicted by many extensions of 
the Standard Model. New color octets can be produced in abundance at the LHC and 
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their production rate is rather model independent. We considered two benchmark models 
as predictive realisations, a simplified sgluon model and a vectorlike confinement scenario 
that features an additional vector resonance in the color adjoint. 

While the QCD backgrounds seem overwhelming at first sight, the use of an equal- 
mass constraint for the two reconstructed dijet resonances and the angular distribution of 
the production of the scalars help to reduce the backgrounds to a manageable level. Both 
models are clearly detectable at the LHC for a broad range of masses. When searching for 
scalar octets with msg > 300 GeV such that m5g/2 is well above the trigger threshold the 
analysis cuts might be adjusted to further optimise the signal yield and increase statistics. 
Comparing the two benchmark models for different masses at low mass only the cross 
section provides some differentiation. If the scalars have higher mass, the reconstruction 
of the four-jet invariant mass might shed light on the presence of an additional s-channel 
resonance decaying into a pair of scalars, however an integrated luminosity bigger than 
lfb~^ is needed. 
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A Cross section ratio 

The ratio of the cross sections for the pair production of hyper-pions and sgluons is given 
in Table 2. The calculation was performed at leading order using the CTEQ6L1 PDFs [46]. 
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Table 2. Ratio of the hyper-pion cross section to the sgluon cross section. 
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